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1 INTRODUCTION 



Recently there has been a considerable amount of progress in understanding 
the physics of two-dimensional quantum gravity coupled to c < 1 matter. 
This was brought by the discovery of the double scaling limit in matrix models 
[Q]: The perturbative expansion of the matrix model provides a discretization 
of the two-dimensional surface, classified according to its genus, which in the 
matrix model is a given order in (A^ being the range of the matrix). By 
tuning the couplings in the matrix model in such a way that the perturbative 
series diverges, one can approach the continuum limit of the two-dimensional 
surface. In this way it is possible to get a sum over all topologies in the first 
quantized string, using as a discretization of the two-dimensional world sheet 
a dense sum of Feynman graphs of an associated matrix model. In certain 
cases, this sum can be calculated exactly, thus providing us with a definition 
of the non-perturbative string. 

Critical exponents have been computed using this method. Furthermore 
it has been shown that a Virasoro algebra naturally arises in these matrix 
models 0]. This algebra is a realization of the Schwinger-Dyson equations in 
the U{N) invariant set of operators of the matrix model. Using this algebra, 
the equivalence of zero-dimensional matrix models to topological 2d gravity 
has been proved P]. Very recently, an application of these ideas to compute 
intersection indices on moduli spaces of Riemann surfaces has been made 

It is clear that further progress in understanding the large- limit will 
provide new grounds to explore the physics of two-dimensional quantum 
gravity and non-critical strings. 

Three years ago we proposed an alternative approach to the large- iV limit 
based upon a hidden BRST symmetry which is present in any quantum field 
theory 0. In particular we were able to prove a theorem that permits the 
solution of the model in leading order of large- in a very simple way . An 
important application of these ideas was made to the two-Hermitian-matrix 
model [^, which corresponds in the string picture to the spherical topology 
of the colored string with two colors 0. In such a case we were able to find 
a set of integral equations that determine completely all correlations of the 
traces of powers of the two matrices, in contrast with previous results where 
only the partition function was found |T^. Very recently, a algebra has 
emerged from this model [jlT] . 

In this paper we study the solutions of these equations in a systematic 
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way. We show that the system of integral equations can be reduced to a set 
of algebraic equations depending on some arbitrary constants. The arbitrary 
constants are fixed by the requirement of having either one or more cuts of 
an associated analytic function. 

Since our main objective is to test the hidden BRST method, we discuss 
in detail several simple examples, the simplest non-trivial case corresponding 
to a cubic potential. We find that the solution of the problem is given by 



the same cubic equation as in ref. [11]. To show explicitly how the arbitrary 



constants are determined, we discuss the cubic potential with one free matrix. 
As a final example we show how to solve the two-matrix model with a quartic 
potential. In this case the system of integral equations reduces to a quartic 
algebraic equation. 

This paper is organized as follows: In section 2 we review the deriva- 
tion of the Schwinger- Dyson equations using the hidden BRST method. As 
an illustration, in section 3 we present a construction of the BRST charge 
encoding the Virasoro constraints of the zero-dimensional, one-hermitian- 
matrix model. In section 4 we derive the integral equations that describe 
the leading order in large- of the Two-Hermitian-matrix model. Section 5 
show how to reduce the integral equations to an algebraic system. Section 7 
contains the examples, and in section 8 we draw some conclusions. 



2 The Hidden BRST method 

In this section we review the BRST-invariant derivation of Schwinger-Dyson 
equations of ref. 0. Usually we study separately the invariances of the clas- 
sical action and the path integral measure. In the hidden BRST method we 
prefer to write a BRST-invariant action that reflects the symmetries of both 
the classical action and the functional measure. To illustrate the point let us 
consider the following functional integral describing the quantum properties 
of a bosonic field (A more general case is considered in ) : 



Z = l [d0]ea;p[-5(0)] (1) 

Here is a scalar field defined in the space-time point x. The func- 
tional measure [d(j)] is invariant under field translations: 
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+ e(x) (2) 

We want to build an action that is invariant under (^. In order to do 

this we introduce a collective field B{x) and write Z (up to a multiplicative 
constant) as follows: 

Z = J[d(f)][dB]exp[-S{(f) - B)] (3) 

S'((/) — 5) is invariant under (|^) provided we transform B at the same 
time by: 



B{x) B{x) +e{x) (4) 

The gauge invariant functions are F[(j) — B] for any F. In the gauge 
B{x) = we get back the original path integral 

In order to quantize the gauge invariant theory that we have just intro- 
duced we use the BRST method[^. The BRST and anti-BRST transforma- 
tions are: 



50 


= ^ 


50 


= ^ 


SB 


= ij 


6B 


= ij 


dip 


= 


6ip 


= -lb 


5i) 


= ib 


6tlj 


= 


6b 


= 


6b-- 


= 



(5) 

To reach the gauge B{x) = 0, we add the following term to the action: 

66[B'^] =6[B^/j] =ibB-i>il) (6) 
If we integrate over the fields b and B the action reduces to: 

^[0] + / [rfx]^^ (7) 
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Using the equation of motion for b and B the BRST and anti-BRST 
symmetry reduces to: 



6S 



6(l){x) 

= f = ee + ee =0 (8) 

We can easily check that(|D is indeed a symmetry of (0). Associated to 
this symmetry is a set of Ward identities, for instance, for any function F: 

< > =0 (9) 



I.e. 



f 6F - 6S 



Computing the (trivial) average with respect to the fermionic variables we 
get: 

The last identity encodes the most general Schwinger- Dyson equation of the 
model. 

Notice that the symmetry @ conmutes with any symmetry of the classi- 
cal action of the model. In particular, if we study a theory which is invariant 
under a large group, the invariance (P) will be respected order by order in 
the expansion. We will use this important fact to derive a BRST for- 
mulation of the Virasoro constraints satisfied by the zero dimensional matrix 
model in the next section. 
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3 Q CHARGE OF VIRASORO CONSTRAINTS 



By now it is well known that the partition function of the zero dimensional 
hermitian matrix model is annihilated by a set of operators that satisfy a 
Virasoro algebra. The existence of this algebra is fundamental to the proof 
of equivalence of two dimensional topological gravity and matrix models 
and to the connection between matrix models and intersection indices in 
modular spaces of Riemann surfaces IQ. 

The Virasoro algebra we are discussing is a direct consequence of the 
Schwinger-Dyson equations of the matrix model. So we expect that the 
BRST-invariant derivation of these equations we presented in the previous 
chapter will produce a ERST extension of the Virasoro algebra. This is 
indeed the case as we are going to explain below. See also (fl^ and||T3||): 

Let us consider the generating function of U{N) invariants correlations 
functions of the zero-dimensional-hermitian matrix model: 

/n=oo 
[dM]exp[- JnTrM^] (12) 
n=0 

Particular critical potentials are obtained by expanding around different 
points in j-space. 

Using the procedure of the last section we form the BRST-invariant ex- 
tension of Z[j]: 

/oo 
[dM] [rfV] [d'ilj]exp[- ^ {jnTrM'^+r]nTrtljM''+f]nTrijM'')]exp[-T 

n=0 

(13) 

Tjn and Tjn are anticonmuting sources. The action Tripip and the integration 
measure are invariant under a nilpotent ERST transformation: 

5M = ^ 5i) = Q # = (14) 

The Ward identity that Z[j,ri,ff\ satisfies is: 

< inriTr(M"- V) -Y^nY. Tr(M"7/'M"~^-™V') >= (15) 



n m 



Replacing 



_^ _dexp\-Tr[^ 
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and integrating by parts over ip gives: 

with the BRST charge Q given by: 

°o n-m S 

^ Of]n+m 



n=-l 



n,m=— 1 



The Virasoro generators L„ are given by: 

d d °° , 9 

n=0 ^3n+m CJn „=o OJn+m 

Notice that due to the identity Trl = we must replace 



-N 



(17) 



(18) 



(19) 



These are the same L„ found in W^. The nilpotency of Q follows because 
the Ln satisfy a Virasoro algebra: 

[Ln,Lj^ = {n-m)Ln+m ri,m = -1,0, 1, . . . (20) 
The same procedure can be used to get a nilpotent BRST charge for the 



Virasoro constraints of the large N vector models discussed in [15 



4 THE INTEGRAL EQUATIONS 

In this section we review the derivation of the system of integral equations 
that describe the large- limit of the Two-Hermitian-matrix model in zero 
dimension 0. 

The action of the model is: 

^ = Tr [Vi(Mi) + V2(M2) - CM1M2] , (21) 

where Ma are N x N Hermitian matrices and Vi(Mj) are arbitrary potentials. 
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We know that S is invariant under Ma UMaU'^UeU{N). Following 
section 2 we form the BRST-invariant extension of S: 

5 = 5 + Tr[V'i^i+^2^2], (22) 

where ipa and i/ja are Grassman valued N x N matrices; S is invariant under 
the following ERST transformation: 

6Ma = ei^a + ei^a (23) 

5tlJa = e6S/6Ma (24) 

6iJa = -e6S/6Ma (25) 

= f = ee + ee = 0. (26) 

As explained in section 2, the Ward identities corresponding to this ERST 
symmetry are the Schwinger-Dyson equations of the model. Moreover we 
can prove the following important proposition, ||^, ^: 

The U{N) invariants of the quantum field are BRST invariants 
in leading order of large-N. 

The proof of this proposition relies heavily on the existence of the exact 
symmetry of both the measure and the action provided by (p6l). 

In what follows, we are going to show how to use this proposition to solve 
the Two-matrix model in leading order of large N. 

Let us write: 

{Mah,=j:<iT:h, (27) 

a 

where are the eigenvalues of the matrix Ma and the corresponding 
projectors satisfying: 

T:T^^ = 5^pT^ E.T- = 1 TrT: = l (28) 

then, since are ERST invariants , according to the proposition stated 
above, we get: 

6{Ma)1^ = E«)'5(Tr)., with k=l,2... (29) 

a 

It follows that the ERST variation of the projectors is: 

ST^ = V " ° + V ° (30) 
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From the proposition, it follows that, in leading order of 1/N , 

STrlT^'^T^^ ■ ■ ■ r°"T^"Va] = (31) 
After computing the fermionic average, we obtain: 

n 1 AT' 

-'ai/3i...afc_ift_i-'7/3fcafc+ift+i...an/3„ _ "r ^ _ 1,<J^; 

We have introduced the definitions: 



A.(.) ^ ^ (33) 

= TrT,-^Tt . . . " (34) 

These relations provide a whole set of identities that determine completely 
the U{N) invariants correlations of the model. Notice that the exact BRST 
symmetry of the model was very helpful in the derivation of the identities. 
Although a complete set of identities follows from the previous equation, in 
this paper we will discuss in detail the case n — 1. 

Taking the continuum version, N ^ oo, of the case n = l,we get: 

Ai{x) -2P J dx'u{x')/{x - x') = c J dyv{y)yF{x,y) (35) 

A^iy) -2P j dy'v{y')/{y - y') = cj dxu{x)xF{x, y) (36) 

P j dx'u{x')F{x\y) / {x — x') 

= [Ai(x) -cy-P j dx'u{x')/{x - x')]F{d-^^) 

P j dy'v{y')F{x,y')/{y-y') 

= [A2(y) -cx-pj dy'v{y')/{y - y')]F{am) 

J dy'v{y')A2{y')F{x,y') = cx (39) 



J dx'u{x')Ai{x')F{x',y) = cy (40) 

(41) 

plus the normalization conditions: 

J dxu{x) = J dxu{x)F(x,y) = 1 (42) 
j dyv{y) = J dyv{y)F{x,y) = 1. (43) 

and P stands for the principal value of the integral, u (v) is the density of 
eigenvalues of Mi(M2) and F{x, y) is the continuum limit of TrT^T2 . Actu- 
ally ( p7D and (|38| ) plus the normalization conditions (0), and (^) imply the 
other equations, as can be easily seen by multiplying by f (y) and integrating 
over y in (^) to get (|35| ) . 

5 SOLUTION OF THE SYSTEM OF INTE- 
GRAL EQUATIONS 

In this section we will show how to solve the integral system of last section. 



Our strategy will follow closely the standard procedure ||T6[. We will explain 
how the system of integral equations reduces to an algebraic equation satisfied 
by an associated analytic function. The different phases of the model are 
determined by the different number of cuts the analytic function may have. 
Let us introduce the following functions: 



U{x,y) = f^' dx'u{x')F{x',y)/{x-x') (44) 

Jd2 

H{x) = f^' dx'u{x')/{x-x'). (45) 

These functions are analytic everywhere in the complex x plane cut along 
the real axis between di and c?2. Across the cut we have: 



(f/+-f/_)/2 = -mu{x)F{x,y) (46) 
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(f/+ + f/_)/2 = P J dx'u{x')F{x',y)/{x-x') (47) 
{H+-H_)/2 = -tTTu{x) (48) 
{H+ + H^)/2 = P J dx'u{x')/{x - x') (49) 

By the index +(— ) we indicate the value of the function in the up- 
per(lower) part of the cut. Moreover, for large x we have: 

U{x,y)r^l/x (50) 
H{x) ~ 1/x. (51) 

According to eq. ( |57| ) we have: 

jj ijj ^ Ai(x) - cy - P J dx'u{x')/{x - x') + iTcujx) 
Ai(x) — cy — P J dx'u{x')/ {x — x') — mu{x) 

But U (x, y) is analytic in the same region of the x complex plane than 
H{x) is; therefore we must have: 

U{x,y)= ,. , /^'''^^ ^, (53) 
[Ai{x) - cy- H{x)\ 

with X{x,y) an integral function of x which is a polynomial in x if Ai is a 
polynomial also; A shall be chosen to get the right asymptotic behavior for 
large x. That is: 

U{x,y) ~ 1/x, for large X (54) 
Plugging ([53|) in (^9]), we obtain 



In a similar way we can introduce analytic functions in the variable y: 



ei 



V{x,y) = / dy'v{y')F{x,y')/{y-y') (56) 

J £2 

/(y) = rdy'v{y')/{y-y'), (57) 



62 
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These functions are analytic everywhere in the complex y plane cut along 
the real axis between ei and 62- Across the cut we have: 



(V+-r_)/2 = -tTxv{y)F{x,y) (58) 

{y+ + V^)/2 = pJdy'v{y')F{x,y')/{y-y') (59) 

(/+-/_)/2 = -t7Tv{y) (60) 

(/+ + /_)/2 = Pjdy'v{y')/{y-y') (61) 

By the index +(— ) we indicate the value of the function in the up- 
per(lower) part of the cut. Moreover, for large y we have: 



V{x,y)r^l/y (62) 

Hy) ~ yy- (63) 

In order to satisfy (^) we must have: 

y ^ ^2{y) -cx-PJ dy'vjy')/ {y - y') + mv{y) 
^ A2{y) - cx - P J dy'v{y')/{y - y') - i7rv{y)' 

Therefore we also have: 

[A2{y) - ex- I[y)\ 

with fi{x,y) being an integral function in y chosen to match the asymptotic 
behavior of V for large y. We get 

F{,x,y) = — — — -. (66) 

[(A2(y) - ex- /+)(A2 -ex- /_)] 

Of course, (|55D and (|66D must coincide for x,y in the support of u{x) 
and v{y). This is a strong requirement, which permits us to find H and /, as 
will be clear after we work out some examples. The main conclusion of this 
section is that the complicated integral system of the last section has been 
reduced to a simpler equality between two functional forms for F{x,y). As 
far as I know this procedure to decouple an integral system was not known 
before. 
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6 EXAMPLES 



In this section we want to discuss some simple models to illustrate the method 
of section 5. We will see that for given Ai and A2, the consistency require- 
ment for F{x,y) will imply a set of equations for the values on the cut of 
the analytic functions we have introduced in the last chapter depending on 
some arbitrary constants. For polynomial A's, this set of equations can be 
reduced to algebraic equations by using the Cauchy theorem of residues. The 
arbitrary constants are the order parameters of the different phases of the 
model. They are determined by the number of disconnected cuts we allow in 
the analytic functions. 

The simplest situation we want to discuss is the "free" case. Namely: 



Ai =x (67) 
A2 =y. (68) 

Actually, we can compute H and / directly by replacing (15) and (16) 
into (^). By the usual method we find: 



u{x) = (1 -c2)Vrf2 -a;V27r (69) 

d = (70) 

v{y) = {l-c')^d^-yy2n. (71) 

Therefore: 



U(x,y) = == (72) 

^ [x{l + c2)/2 -cy + {l- c^)Vx^ - dy2n] ^ ' 

" [y{l + c2)/2 - cx + (1 - c^)^y^^/2Ti] ' ^^^^ 

The square root takes its principal value with a cut along the real axis 
between —d and d. It follows that: 

" [C2(x2 + y2) + l_c2-c(l+c2)xi/]' ^^^^ 
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It is easy to check that ( |7^ ) gives the right perturbative expansion in c. 
Let us examine a more general case: 



Ai(x) =x + gix'^ (75) 

^2{y) =y + g2y^ (76) 

According to ( [5B| ) and (^3]) we must have: 

X + gix^ — cy — n[x) 

V(.,y)= « + (78) 

The comparison of the F{x,y) coming from U{x,y) and from V{x,y) 
gives: 



/iic^ — c(72(2Ai — — HJ) = a2X^ + aix + (79) 
^72(Ai - iJ+)(Ai - H^) - c/ii(2Ai - - if_) 

= 622;^ + bix + 60 (80) 

/ii(Ai - i7+)(Ai - i^_) = 5(ic^x^ + /i2X^ + /iix -(M^ 

and 

Aic2-c(7i(2A2 -/+-/-) = a2l/' + + /i2 (82) 
^2(A2 - /+)(A2 - /-) - cAi(2A2 - /+ - /-) 

= aiy^ + biy + hi (83) 

Ai(A2-/+)(A2-/-) = g2cV + +aoy^ + boy -m 

for certain constants ai,bi,hi. 

From equations (11) and (12), we get: 

/ii = g2[A,-{H+-H_)]/c+l, (85) 
Ai = (7i[A2-(/+-/-)]/c+l, (86) 
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from which we can fix several of the arbitrary constants: 



a2 = -gig2C (87) 

ai = -92C (88) 

Co = (89) 

h2 = -9ic (90) 

h2 = (91) 



The remaining system of equations is: 



Hi + H+H_ + Hl- {H+ + /f_)(2Ai + 0/92) + 

A? + 2Aic/^2 + fi = (92) 
((?2(Ai - (i7+ + /f_))/c+l)[Ai -/f+][Ai = ^icV + cV + /iix -(0^^ 

where we have defined 

B = {-gic^x^ + hx + bo)/g2 (94) 



To solve this system, multiply eq.(p^) by (H^ — HJ). A remarkable property 
of the resulting expression is that it depends solely on the discontinuities of 
powers of H{x) along the cut. Using this property, we can reduce the system 
(|92D and (p3| ) to a cubic equation for H{x). In fact, consider the integral: 



1 r , Hixf 
dx- 



27rz J X — y 

along a curve that encloses the cut of H{x). In the region of the complex x 
plane that does not contain the cut, the integral equals H{yY. On the cut, 
we can use ( p2| ) to express H+{x)^ — H_{x)^ as a function of H^{x)^ — H^{xY 
and H+{x) — if_(x). Then the remaining integrals can be computed in the 
region that do not contain the cut to get the following equation for H: 



H{yf - (2Ai + c/g2)H{yf + F{y)H{y) - F'{y) - F"{y){< x > -y)/2- 
F"'{y){< x'^ > -2 < X > y + y^)/6 - F^^\y){< > -3 < > ?/ + 

3 < x > ?/2 _ y^)/2A + 2gi = 0(95) 
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for certain additional constants < >. We have introduced the notation: 



F{y) = A? + 2cAi/(72 + B (96) 

The condition ( [^ ) has been exphcitly used to derive (P^D, which coincides 
with equation (4) of ref. [|ll] . 



Let us review what we have obtained up to here. The first step in the 
solution of the set of integral equations of section 4 was to express f/(a;, y) 
{V{x, y)) in terms of H{x){I{y) plus certain arbitrary functions of x{y). The 
next step was to impose the equality 

U+{x,y) -U-{x,y) ^ V+{x,y) - V_{x,y) ^ x /pyx 

H+{x)-H4x) I+{x)-l4x) ^ '^^ ^ ^ 

In this way, we found conditions that relate these arbitrary functions. Re- 
markable enough, these conditions can be solved completely, reducing the 
arbitrariness to a set of constants. This is the main result of this paper. It is 
clear that the same procedure will work for an arbitrary potential. 

In what follows, we will sketch that solution of the cubic equation which is 
relevant for the phase of the Two-matrix model containing the perturbative 
solution. That is, we will explore the conditions under which H{x) has only 
one cut on the real axis. This requirement will fix all the arbitrary constants. 
In principle, we could find all the phases of the model allowing more cuts, 
but we will not do this here. 

6.1 Analysis of the cubic equation 

The cubic equation + 02-^^ + aiH + ao = has the following solutions: 



Hn = WnSl + wls2 - 02/3 (98) 

S, = [r + (r2 + g3)l/2]l/3 ^gg) 

S1S2 = -q (100) 

q = aj3-{a2/3f (101) 

r = (aia2 - 3ao)/6 - (02/3)^ (102) 



where Wn is the n-th root of unity. 
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The equations ( p2D and ( ]93| ) are satisfied if: 

q = 5/3 + (Ai - c/(72)V9 (103) 

r = c{gic^x^ + c^x^ + /iix + ho)/{2g2) + g(Ai - c/^2)/2 + 
(Ai - c/(72)V54 (104) 

bi = gig2 - cic'' + 1) (105) 

Tlie arbitrary constants will be determined by demanding that the ana- 
lytic continuation of the solution of the cubic that behaves as 1/x for large 
X must have only two real branch points. 

In our case this requirement is satisfied if r"^ + q^, which is a polynomial 
in X of degree 9, can be written as follows: 

+ = —A{x^ + p2x'^ + pix + Po)'^{x^ + S2X^ + SiX + Sq) (106) 

for some constants pj, Si and A. 

Actually we can have different solutions of the problem by permitting a 
larger number of branch points , but we will not consider them here. In 



the next subsection we study a simpler limiting case of the last equation. 
6.2 One of the matrices is freer^^i = 

In this section we consider the limit gi = The equation for I{y) reduces 
to a quadratic equation, which can be solved easily. The one cut solution 
implies the following values for the arbitrary constants: 

/ii = (107) 
bbo = 1 (108) 

and ho satisfies the following equation: 



64glhl + hl{96c'gl - 96gl - (c^ - 1)^) + 
hoiSOc'gl - mc^gl + W2 - (c' - 1)') + 

21gt-c''gl + ?,c^gl-?,^gl + gl = 0. (109) 

Now we want to check the consistency of the whole procedure by finding 
hi, bbo, ho using (|106|) . For = 0, + is a polynomial of degree 5; in 
order to get just two branch points, we must thus have: 



17 



+ = -A{x + po^ix^ + S2X^ + six + sq). 
We find the following equations to determine po and ho: 

ho = -{c^g2 + 0^92 - 5cg2^ - Gc^gs^ + c^Po - 2c^Po + c^Po - 

^c^g2^Po - '^c^92^Po - ^'^c^g2^Po + g2'^Po - ^c^g2Po'^ - 

3c''g2Po^ + l2c^g2Po^ - Gc'g2Po^ + 3cg2^po^ + 30c^g2^po^ + 
2c^g2''po' + lQc''g2'po' - IQc^gz'po'' - 10c=^^/p/)/ 
(3c=^(?2 + 9c^^2 + 9c^/ + Qc^g2^Po - 3Qc'g2^po - 6cg2'po^), (HO) 



= c^g2 + ^c^g2 + 3c^^2 + + (c^ + - 5c^ + 3c^ - 12c^^/) ])o + 

The last two equations for po have a common root if ho satisfies ( |109D . 

As usual the asymptotic expansion of the functions H{x) and I{y) around 
oo gives Green's functions of the eigenvalues of the corresponding matrices. 
For instance, we get: 

c{-l + c^ + ho) 
< X > = — 

g2 

^ 2^ c'-2c' + c^ + g2'-c'ho + c'ho 

< X > = 5 (1121 

gi 

When g2 7^ 0, the system of equations to determine the constants in 
( |106D turns out to be too complicated to be solved analytically, so we will 
not discuss it here. 

Instead we will mention how these results generalize to a quartic potential. 



6.3 Quartic Potential 

For the sake of completeness, we will sketch the solution of the two-matrix 
model with a quartic potential, which is defined by: 

Ai=ax + gix^ (113) 
A2 = by + g2y^ (114) 
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So, we must have: 



U{x,y) 



(115) 



gix 



+ ax — cy — H{x) 
92y^ + fJ-iy + 1^0 



(116) 



92y 



+ by — cx — I{y) 



Again, imposing the equahty of F{x, y) coming from the two previous 
functions, we get a system of algebraic equations to determine the arbitrary 
functions Aq, Ai, /io, /ii-One of these equations is: 



g^iHX - Ht) - (3(72Ai + c){Hl - Hi) + {?,g^^l + + 3cAi)(i/^ - Hi) - 




This system can be solved using the Cauchy theorem of residues as was 
used to derive (|95|) . We get a quartic equation for H{x) and a similar one 
for I{y). 

Again the complicated set of integral equations has been reduced to an 
algebraic equation, depending on a set of arbitrary constants. To completely 
determine the arbitrary constants we may demand that H{x) has only one cut 
in the x complex plane. Once this choice has been made all U (A^)-invariant 
correlations of the two matrices which involve u{x),v{y) and F{x,y) can be 
computed. 



7 CONCLUSIONS AND OPEN PROBLEMS 



We have been able to solve the system of integral equations that describe the 
large- limit of the two-matrix model( a particular subset of the Schwinger- 
Dyson equations) by reducing it to an algebraic equation satisfied by an 
associated analytic function which depends on some arbitrary constants. In 
particular, we have verified that the solution found in ||ll| is obtained. In 
certain simple cases, we have been able to find these constants explicitly by 




Ai{g2Al + 2B + 3cAi)(i/+ - H^) = Dc{H+ - HJ)m) 



with 



B = b2x'^ + bix + bo 



and 



D 



gicbx^ + d2x'^ + dix + cLq 
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choosing the analytic function with only one cut (the perturbative phase of 
the model). In more general situations, finding the constants for the one cut 
solution analytically is very complicated. 

Our solution is appropriate to explore all the different phases of the model 
which are characterized by the set of arbitrary constants. The constants are 
no longer arbitrary if we fix the number of branch points of the associated 
analytic function. We see this result as a strong test of the hidden BRST 
method. 

It is an interesting open problem to find the other solutions of (32), which 
will enable us to compute more general correlations functions of the form 
TrM'l^Ml^^M'l^Ml^^ .... 
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